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' We use the LUQAS sample of 27 high resolution, high signal-to-noise QSO absorption spectra 

, (Kim et al. 2004) and the results from Croft et al. (2002) together with a suite of high- 

""■■-^ ■ resolution hydro-dynamical simulations run with the GADGET-II code, to infer the linear 

' dark matter power spectrum on scales of ~ 0.3 — 30 Mpc at z — 2.125 and at z = 2.72. 
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The Lya forest in QSO absorption spectra has been used to probe dark matter power spectrum 
a number of authors (e.g. Croft et al. Gnedin & Hamilton^, Viel, Haehnelt & Springel 
and references therein). The Lya forest is sensitive to fluctuations in the DM density on 
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scales of a few Mpc. It can thus be used to investigate the possible cut-off of the DM spectrum 
which is expected if the dark matter were warm dark matter, it can give constraints on the 
matter fraction in neutrinos and allows to investigate the gravitational growth of structure and 
possibly the redshift evolution of dark energy (Viel et al. EB, Mandelbaum et al. ^ Lidz et 
al. 13). Verde et al. ^ combined the Lya forest results of Croft et al. with the CMB results of 
WMAP and concluded that there is evidence for a running spectral index and for a tilt in the 
primordial power spectrum. However, a recent analysis by Seljak et al. El argued that the errors 
are larger and the effective optical depth is smaller than that assumed by Verde et al. Seljak et 
al. concluded that there is no evidence for a running spectral index or for a tilt of the primordial 
power spectrum. We briefly discuss here the flux power spectrum of the LUQAS sample (Kim et 
a set of high resolution QSO absorption spectra, and compare it to the recently published 
flux power spectrum obtained for a large sample of low-resolution SDSS spectra (McDonald et 
al. ™). We will further discuss the linear dark matter power spectrum which Viel et al. 
obtained from the LUQAS flux power spectrum using the effective bias method proposed by 
Croft et al.*i. Implications for the rms fluctuation amplitude as and the primordial index n are 
also discussed. 




Figure 1: Comparison between the flux power spectrum of the LUQAS sample (continuous line with error bars) 
and the recent results for a large sample of SDSS spectra by McDonald et al. 2004 (dashed line). Note that the 
SDSS data refer to their table 3 including the background contribution. 



2 The LUQAS sample 

The LUQAS sample (Large Sample of UVES Quasar Absorption Spectra) consists of spectra of 
27 QSOs taken with the Ultra- Violet Echelle Spectrograph (UVES) on VLT. Most of the spectra 
have been taken as part of the Large ESO Observing programme UVESLP (P.I.: J. Bergeron). 
The median redshift of the sample is z = 2.25 and the total redshift path isAz = 13.75. The 
typical signal-to-noise ratio is ~ 50 and the pixel size is 0.05 A. For a more detailed description of 
the sample and the data reduction we refer to Kim et al. In Figure 1 we compare the ID flux 
power spectrum of four subsamples of LUQAS at different redshifts (continuous curves) with 
the flux power spectrum obtained by McDonald et al. ^ from a large sample of low-resolution 
SDSS spectra (dashed curves). 

There is reasonable good agreement between the two estimates of the flux power spectrum 
over a wide range of wave-numbers, in all four redshift bins. Note that due to the significantly 
higher resolution the LUQAS sample can constrain the matter power spectrum all the way to 
the thermal cutoff (k > 0.01 s/km). At these small scales the estimate of the SDSS sample is 
affected by its low resolution. The statistical errors of the SDSS sample are significantly smaller 
but the systematic uncertainties are larger than the statistical errors (Kim et al. McD onald 
et al.Sl). The flux bispectrum of the LUQAS sample has been presented in Viel et al. 02. Note 
that the large discrepancy between the two samples claimed in McDonald et al. was mainly due 
to an an erroneous shift of the flux power spectrum by half a bin-size in log k in Kim et al. 

"Tables of the flux power spectrum are available at: www.ast.cam.ac.uk/~rtnigm/luqas.html. Please note that 
the k- values of tables 3 and 5 of Kim et al.^had been erroneously shifted by half a bin-size in log k (erratum in 
press). Corrected tables can also be found in astro-ph/0308103 
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Figure 2: The recovered linear dark matter power spectrum at two different redshifts. At z = 2.f25 using the 
LUQAS sample (left panel) and at z — 2.72 using the Croft et al. (2002) sample (right panel). The filled 
triangles have been obtained with an effective optical depth derived from high resolution spectra and from hydro- 
simulations with To = fO 4 ' 15 K and 7 = f.6. In the right panel we overplot the original result from Croft et 
al. (empty triangles) and a re-analysis of their data with the effective optical depth they assumed (diamonds). 
Continuous black line shows a model which is a good fit to our data. The dashed line represents a model with 
a running spectral index with best fitting values suggested by the WMAP experiment (Table fO, Spergel et al. 
2003). The yellow rectangles indicate the range of wavenumbers used in our analysis. 



(2004, erratum in press). 

3 Inferring the linear dark matter power spectrum 

The method we use to infer the linear dark matter power spectrum has been proposed by Croft 
et al. - (but see Gnedin & Hamilton^). It uses numerical simulations to calibrate the relation 
between flux power spectrum and matter power spectrum. It then assumes that the flux power 
spectrum Pp(A;) at a given wavenumber k depends linearly on the linear real space matter power 
spectrum P{k) at the same wavenumber and that both can be related by a simple bias function 
b(k) as P F (k) = b 2 (k)P(k). 

We have run a suite of hydro-dynamical simulations with varying cosmological parameters, 
particle numbers, resolution, boxsize and thermal histories using a new version of the parallel 
TreeSPH code GADGET, in order to explore the different systematics and statistical errors. 

In Figure 2 we show the results of our analysis in terms of A 2 (k) = P(k) /c 3 /(27r 2 ) for the 
LUQAS sample (left panel) at z = 2.125 and for the Croft et al. sample (right panel) at z = 2.72. 
The filled triangles represent the recovered linear dark matter power spectrum using values of 
r c fj suggested by high-resolution spectra. 

The dashed curves in Figure 2 represent the best fitting running spectral index model of 
Spergel et al. El which is a fit to WMAP, CBI, ACBAR, 2dFGRS and the Lyman-a forest data, 
while the continuous line is a good fit to our data points. 

In order to minimise uncertainties due to continuum fitting and metal lines, and to avoid 
dealing with the problematic thermal cut-off at small scales, we have only used the range of 
wavenumbers 0.003 < k (s/km) < 0.03 for our quantitative analysis, which corresponds roughly 
to scales 0.3-30 Mpc. 



The main results can be summarized as follows. With the same assumptions for effective 
optical depth, density-temperature relation, and cosmology, our inferred linear matter power 
spectrum (empty diamonds right panel) agrees very well with that inferred by Croft et al. ^ 
(empty triangles right panel). 

We confirm previous results that the inferred rms amplitude of density fluctuations depends 
strongly on the assumed r c ff. It increases by 20% if we assume an optical depth of T e s = 0.305 
a value indicated form high-resolution absorption spectra. We find, however, a dependence on 
r^ff which is weaker than that of Croft et al. - and stronger than that of Gnedin & Hamilton 
121 . The decrease of the amplitude of the flux power spectrum between z = 2.75 and z = 2.125 
is consistent with that expected due to the decrease of T e s and the increase of the amplitude of 
matter power spectrum due to gravitational growth. 

Our estimate of the systematic uncertainty (error bars in the bottom part of the two panels 
in Figure 2) of the rms fluctuation amplitude of the density (~ 14.5%) is a factor 3.5 larger than 
our estimate of the statistical error (~ 4%). The systematic uncertainty is dominated by the 
uncertainty in the mean effective optical depth and by the uncertainties between the numerical 
simulations of different authors. Reducing the overall errors will thus mainly rely on a better 
understanding of a range of systematic uncertainties. 

By combining the constraint on the amplitude of the DM power spectrum on large scale from 
COBE (assuming that there is no contribution from tensor fluctuations) with the constraint from 
the high-resolution Lyman-a forest data on small scales we obtain n = 1.01 (O 0m /i 2 / - 135 ) -0 ' 35 ± 
0.02(statistical) ± 0.06(systematic) for the spectral index. The corresponding rms fluctuation 
amplitude is, a% = 0.93 ± 0.03(statistical) ± 0.09(systematic) . 

Thus, for values of the mean optical depth favoured by high-resolution spectra, the inferred 
linear power spectrum is consistent with a ACDM model with a scale- free (n = 1) primordial 
power spectrum. 
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